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ABSTRACT 


It is well known that conventional design practice 
of RGG members often results in overly conservative designs. 
It is also known that the degree of conservatism can not 
readily be estimated by current methods. The statistical 
nature of design variables is usually ignored in conventional 
practice. The growing need for answers to such questions - 
How safe is design ? How can designs be made to specified 
levels of adequacy (reliability)? - leads directly to 
probabilistic methods. The present work is a step towards 
this direction. 

This work can be broadly classified into two parts. 

(i) Reliability based design of RGG circular slabs with 
fixed as well as simply supported edge conditions. 

(ii) Optimum cost design of RGG circular slabs with a 
reliability constraint. 

The first part of the work consists of the parametric 
study of the reliability based design of under “reinforced 
circular slabs with uniformly distributed load. By changing 
the value of each parameter turn by turn the effect of design 
variables on the area of steel reinforcement is studied. 

The design parameters included in parametric study are depth 

/ 
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of the slab, characteristic compressive strength of concrete, 
characteristic yield strength of steel, reliability parameter 
and coefficient of variation of the design variables. It is 
found that the area of steel is very much sensitive to the 
yield strength of steel and reliability paramter. 'The steel 
area is also sensible to the coefficient of variation. 3y 
grouping the parameters in the nondimensional form, reliability 
analysis is carried out and design curves are presented as per 
Indian Standard Code Practice (IS 456-1978), Two examples are 
taken for illustration of the use of the design charts. 

Finally the optimum cost design of under reinforced 
concrete RCO circular slabs is studied with a reliability 
constraint. The design is considered safe and adequate if 
the expected moment capacity is greater than design moment 
capacity by a certain number of standard deviations. The 
effect of coefficient of variation and change in cost ratio 
on the optimum design cost is also studied. The optimal 
cost design problem is cast as a nonlinear mathematical 
programming problem and is solved by using sequential 
unconstrained minimization technique. 



CHAPTER 1 

INTRODUCTION 


1 . 1 General : 

'The fundamental role of probability theory in safety 

and performance analysis is widely recognized in all branches 

of engineering. Despite more than 20 years of recognition in 

structural engineering literature , the probabilistic methods 

have not yet been explicitly adopted as a basis for any standard 

code of practice. The probability theory provides a more 

accurate engineering representation of reality. Many leading 

(1 ) 

civil engineers in several countries have written about 
the statistical nature of loads and material properties. It 
has been demonstrated that such uncertainty in applied forces 
and in structural resistance implies uncertainty in structural 
performance and that this uncertainty can be analysed rationally 
only with probability theory. 

The conclusion is that, if structural safety is to be 
placed in a position where it can be discussed quantitatively, 
it must be treated probabilistically. Several authors suggest 
or imply that the traditional nominal strengths, nominal loads 
and safety factors (or load factors) be replaced by an allowable 
’probability of failure 1 calculated on the basis of observed 


frequency distribution of strengths and loads. In this way 
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it is claimed that more realistic statistical models would 
produce better engineering results. 

'Reliability' is the capacity of a structural system 
to perform its assigned functions under specified environmental 
conditions for a specified period of time. In otherwdtds, it 
is a measure of the adequacy of a design in its intended environ- 
ment, whenever physical or other factors that must be considered 
are essentially probabilistic. 

1.2 Reasons for Slow Inclusion in Oodes of Practice: 

Although a considerable amount of research effort has 
been devoted to the study of probabilistic approach of design, 
the progress towards its implementation in standard codes of 
practice is rather slow due to the following reasons: 

(i) Por numerous reasons standard codes must remain uncom- 
plicated in form and must have a basis that is straight forward 
to understand. The reliability analysis calculations in the 
literature v do not promise either benefit, Aitnough a new code 
format that the profession believed to have a more rational 
basis and to lead to more consistent safety could presumably 
afford .some minor increase in complexity over present codes. 

(ii ) The assessment and Justification of a numerical value 
for the 'acceptable risk' of failure. 
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(iii) Engineers recognize that present safety factors 
account for many uncertain! ties which are seldom included 
in the observed statistical data. 'These other factors 
include those about which there is incomplete knowledge and 
those for which the engineering analysis is possible but not 
economically justified (e.g. the relationship between labora- 
tary strength and field strength ; local ,long term and non- 
linear stress distributions; actual versus assumed uniform 
load distribution, perhaps even minor engineering miscalcu- 
lations and construction blunders whose guranteed elimination 
would require many hours of tecnnical man power). In short, 
much of the uncertainty in structural performance predictions 
is of a fundamentally non statistical nature. 'This is not to 
say however, that it is nonprobabilistic in nature. 

(iv) On a more esoteric level, it is unreasonable to 
believe that a complete, strictly statistical probabilistic 
analysis of structural safety will even be feasible for 
standard structures. The available data will always be too 
limited to provide wholly reliable estimation of failure 
probabilities of the order of magnitude necessary for adequate 
public safety. The sample sizes required are in multimillions. 
Extrapolation of mathematical models of frequency distributions 
will always be a part of these analyses. 

There are, then sound reasons why the suggested adoption 
of probabilistic safety analysis has not been implemented in 
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standard. structural engineering practice, Nonetheless, the 
fundamental' arguments for a statistical basis remain. It 
still promises to provide a more rational, quantitative 
representation of engineering design. 

1.3 fallacies in Designing by Safety factors: 

The quantity in. use to maintain a proper degree of 
safety in a conventional design is usually referred to as 
the factor of safety. The following are the several other 
definitions in its broader sense, i.. e. ratio of the maximum 
safe load to the normal service load, ratio of mean strength 
to mean load or mean failure governing strength to mean 
failure governing stress. 

There are several important, relatively new. circum- 
stances which have brought the inadequacy o.f, .the factor of 
safety concept to the fore. Some of them may be mentioned as : 

(1 ) The statistical nature of design variables is usually 
ignored in conventional practice. The conventional ' approach 
in design practice may be compared to a kind of worst-case 
analysis. The maxima of loading and the minima of strength 
are treated not only ah representative of design situations, 
but also of simultaneous occurrence. This is the^ basis on 

...which unknown parameters 'are' computed. .Actually, magnitude 

'* * * * 

,.$nd frequency relationships for both load and strength must 
be considered to avoid unrealistic results. If an extremely 
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largo load must act on an extremely low value of strength 
to induce a failure, then the probability of such simultane- 
ous occurrences is very important. Fig. 1.1 picture the 
type of comparison discussed above. 

(2) The safety factor' concept completely ignores the 
facts of variability that result in different reliabilities 
for the same safety factor. The Fig. 1.2 indicates a finite 
probability of failure whose magnitude is a function of the 
degree of overlap of the two distributions. As the overlap 
increases, the shaded area and consequently failure proba- 
bility increase proportionately. Thus there is a rational 
explanation for changes in failure probability which 
accompany changes in the load -strength distribution overlap. 

Three possibilities exist in which a safety factor (SF) 
may be maintained as failure probability varies; 

(i ) Mean load and mean strength may be changed in the 
same proportion with no change in the standard deviations. 

Thus , 

R K-|R 

ST? = — = -1- = K ... (1.1 ) 

L K t I 

where K-j'R and K-jl reflect a shift in mean strength and load 
respectively. 'The Fig. 1.3(a) illustrates the shift for 

> 1 . This demonstrates the limitation of safety factor as 
an indicator of failure incidence and of design safety. 
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(ii) Failure probability varies with the mean values and 
distributions are held constant and the standard deviations 

i 

are varied as shown in Fig. '1.3(b). Since the mean values 
are not changed, the safery factor remains the same, but 
failure probability does change. 

(iii) It is possible to change both mean and standard 
deviation values without changing the safety factor. 

Fhe Table I clearly indicates that even though the 
factor of safety remains same, the probability of failure 
may vary significantly. 

TAB IE I 

SAFETY FACTOR M D RELIABILITY 


Mean 

Yalue 

Standard 

Deviation 

Safety 

Factor 

Reliability 

Strength 
R(MPa ) 

Load 

L(MPa) S 

Strength 

> R (MPa) 

Load 

S^MPa) 

250 

100 

10 

15 

2.5 

0.9166 

250 

100 

80 

30 

2.5 

0.9099 

250 

100 

100 

60 

2.5 

0.8997 

250 

100 

250 

250 

2.5 

0.6628 


Bar (-) over the variable indicates their mean value. 
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(3) 'The day to day increasing complexities and environmen- 
tal actions presenting a much broader spectrum of loading 
(for 3-d structures, concrete structures, marine structures, 
nuclear weapons and space explorations etc. ) necessiate the 
study of statistics of extremes for the rational design of 
such structures. 

1.4 Need for Probabilistic Approach to Reinforced Concrete 
Members : 

The actual strength of a reinforced concrete member 
(RCC) differs from the nominal strength calculated by the 
design engineer due to variations in the material strengths 
and geometry of the member, as well as the variabilities 
inherent in the equations used to compute the member strength. 
Similarly, designers use constant nominal values of loads 
in their calculation of forces, but the actual loads are 
variable. This variability in strength and loading is 
accounted for in one form or another in safety provisions 
of all existing codes. As a first step in the calculation 
of resistance factors, it is necessary to study the variability 
of the strength of reinforced concrete in flexure, shear, 
combined axial force and flexure etc. This requires a 
knowledge of the variability of the parameters that effect 
the strength of RCC members. These are the concrete strength 
in compression and tension, the yield strength of steel and 
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the dimensions of the cross s'ection of the member. In 
practice the compressive strength of concrete varies with 
in the same mix also. Por considering this variability 
in compressive strength of concrete IS code of practice ^2) 
is specifying characteristic strength (i.e. more than 5 
percent of the test results should not fall below the 
specified strength). Even at the time of construction 
also the cross sections of the members may be slightly 
different from section to section because of the construction 
faults . 

In reliability approach all these variabilities are 
taken in to account. Because of these variabilities in 
materials and material strengths of RCC members, it is 
more appropriate to go for a reliability approach. 

1.5 Steps Involved in Probabilistic Analysis and Design: 

The various steps involved in probabilistic analysis 
and design may be presented in a flow chart as shown in 
Pig. 1.4. 

1.6 Scope of the Thesis : 

The major interest of the thesis is to study the 
design aspects of RCC circular slabs with preassigned 
probability of failure. The work can be broadly classified 
into three parts. 
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i) Parametric study of the design formulation for RCO 
circular slabs, finding the area of steel by changing a 
particular design variable while all other variables are 
kept constant. 

ii) Presentation of some design charts for RCO circular 

slabs with reliability parameter as per Indian Standard 

( S2 ) 

Code of Practice for Plain and Reinforced Concrete . 

iii) Finding the minimum cost of RCC slabs with reliability 
as a constraint. 

In the first two chapters, the general aspects of 
reliability based designs are reviewed, fallacies in safety 
factor concepts discussed and literature is reviewed. 

Ihe Chapter 3 contains design formulation for RCC 
circular slab design and some readymade design charts for 
free edge conditions as well as for fixed edge conditions. 

The Chapter 4 deals with some studies on optimum cost 
of the RCC slabs with reliability as a constraint. 

Some design examples, discussion and conclusions are 


presented in Chapter 5 



o - L R L,R 

FIO-1-2 FREQUENCY DISTRIBUTION OF THE LOAD (L) 
AND MATERIAL RESISTANCE (R) INCLUDING 
PROBABILITY OF. FAILURE (Shaded area) 






(b) 


FIG. 1-3 EFFECT ON FAILURE PROBABILITY DUE TO 

CHANGE IN MEAN AND STANDARD DEVIATION 
VALUES 
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Pig. 1.4 Steps Involved in Probabilistic Analysis and Design 
















CHAPTER 2 


LITERATURE REVIEW AUD REII ABILITY THEORIES 
2.1 literature Review: 

Only in recent years a sizeable quantity of literature 
has appeared in the area of reliability based design and 
probability of failure of structures due to random variables 
of loads and strengths. 

Literature on probabilistic distributions for strength 
of concrete beams, safety of structures, risk analysis, pro- 
babilistic design of members and optimum designs with relia- 
bility constraints is reviewed. • 

Among the first to recognize the concepts of an 

associated probability of failure in any design are Ereudenthal 

( p ) ( E ) 

Pugsley v 1 and Prot . Preudenthal analyzed the safety 
factor in engineering structures while Pugs ley pleaded for 
the assessment of probability of failure and the collection 
of statistical data of loads and structural properties. 

(1 ) 

Preudenthal analyzed the safety factor in 
engineering structures in order to establish its magnitude. 

He derived the safety factor from absolute and measurable 
physical properties and phenomenon. External conditions and 
pertinent physical properties and phenomenon were analysed 
under the common aspect of their relation to the magnitude 
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of safety factor. Pugsley^ discussed the method of design 
incorporating the safety assumptions as well, as the efficiency 
concepts. ihe terms like structural efficiency and structural 
safety were defined and a strut problem was considered for 
the classification of concepts. . 

( 4 ) i:: ' 

Freudenthal outlined the procedure for the 
numerical evaluation of factor of safety based upon the 

' 0 OI, " 

concept of unserviceability or probability of failure. 

Preudenthal et.al. J correlated the factor of 
safety with probability of survival or failure. It also 

provides guidance by suggesting and illustrating techniques 

' 0 u;c» :;1 

that may be of considerable value in studying certain phases 
of the safety (reliability) problem. 

Een.jamin^ discussed the superiority of probabilistic 

0 •; • ; ■ ) i..v & !•; y 0 1 ' 

designs over conventional deterministic procedures in respect 
of informational content, modeling of reality, refinement of 
analysis and design, decision making and the resulting 
structures. 

( 7 \ 

Sexsmithf J presented procedures for determining the 
probability distribution of the safety margin and related 
quantities for several types of reinforced concrete members. 
Data were extracted from previously published test results, 
and were separately analyzed for each of three common 
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structural types; beams without compression reinforcement, 
beams with compression reinforcement and tied columns. 

(s ) 

Benjamin and Iind discussed one possible application 
of probabilistic procedures to back-up a deterministic code. 

Subjective probability measures were used to obtain model 

# 

parameters and then some applications of the results were 
presented. 

Shahh outlined reasons for interest in probabilistic 
studies of safety principal considerations in formulating 
probabilistic codes were discussed and a simple strength 
modelling example was introduced. The probability impli- 
cations of ACI 3 18-6 3 ultimate strength design procedures 
were presented. 

Cornell introduced a frame work for structural 
code in which probability was used to enhance realism. He 
treated all the uncertain! ties through standard deviation 
and presented a format of AGI code of practice (AGI 318-63). 

, ( 1 1 ) 

Alien' ' has summarized the results of a probabilistic 
study of ultimate moment and ductility ratio of reinforced 
concrete. The results presented are useful in assessing 
existing design formulas and safety provisions. 

(ip) 

Moses and Stevenson J presented methods for 
incorporating reliability analysis into optimum design 
procedures. Examples were given for multimember elastic 
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truss design and frames designed according to limit design 
theory. The approach adopted was to design for a specified 
probability of failure, in which the failure probability 
whs evaluated from a sequence of numerical integrations. 

Ang v J has summarized some basic probabilistic models and 
he had shown that practical formulation of risk based design 
and evaluation of safety and performance, including repeated 
loading and fatigue can be developed consistently within 
the existing probabilistic models of structural reliability. 
Chandrasekhar ^ ^ has analyzed number of sets of field data - 
and found that the probability of failure of concrete cube 
is about 6.12 for' the strength ranging from 35 MPa and the 
individual variation of strength from the mean value is 
+ 10 percent. h r 

kng^^ identified the requirements and conditions 
necessary for implementing the ^ibhabilistic approach and 
has shown that the explicit consideration of risk in the ... ■ 
evaluation of performance, and formulation of design are now 
practically feasible and desirable. ncr-- .- 

n / 1 f 7 S ' Y r 

Ranganathan and Dayaratnanh _ used Monte Carlo 

' ( ..'I . , '' • ,.;y ; 

Simulation fechniq.ue- to ...evaluate the probability of failure 
of a prestressed’ concrete flanged section. The failure of 
the section at limit state of strength was divided into two 
cases: (i) under-reinforced and (ii) over-reinforced and 
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the probability of failure was evaluated based on the above 
two conditional failure events. 

( 18 ) 

Shah and Sexsmith ' described the general philosophy 
of developing probability based codes. Various probability 
based improvements were suggested for the current ACI Code of 
practice. 

( 2 ) 

Mirza et.al. reviewed the variability of strength 
and stiffness normal -weight structural concrete and 
representative distributions were suggested for use in 
estimating the effect of these variations on the strength 
of reinforced concrete elements. 

(21 ) 

Bllingwood described the results of a study 
where objective was to assist in the development of probabi- 
listic limit state design criteria for use in the United 
States of America. 

( 22 ) 

Mirza and MacGreger documented the assumptions and 
analysis made in the derivation of resistance factors for 
reinforced concrete design compatible with the load factors . 
in the National Building Code of Canada. Statistical estimates 
were made of the variability of the strength of reinforced 

1. 

concrete members based on the available data on the variability 
of' concrete and steel dimensions. 

MacGregor et.al. presented the results of studies 

of the variability in the strengths of reinforced concrete 


U'vVv 
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beams and columns. They also presented the coefficients of 
variation of depth and steel reinforcement of slabs. For 
variabilities of member strengths, they had interpreted 
the results of Monte Carlo Simulations; these inturn were 
based on statistical distributions obtained from literature. 

Roy(^) presented some design charts for reliability 
based design of RGC5 beams. He has taken the coefficients of 
variation for different design variables from test results 
compiled by -Allen on more than 100 beams. Finally optimiza- 
tion of beams with reliability constraint was presented. 

( 25 ) 

Moses and Kinser presented optimum sizing of 
elements for a multielement, multiload conditions structure 
with safety defined interms of allowable probability of 
failure. 'They presented weight minimization and probability 
analysis which are applicable to any elastically designed 
structure and can treat any frequency distribution for each 


loading and element strength. 

Rao^^ discussed the optimum cost design of under- 
reinforced concrete beams by treating all the design parameters 
as random variables and considering the reliability as the 
constraint, numerical 1 examples were presented to illustrate 
the effectiveness of the method. 
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2.2 General Expression for Reliability: 


For deriving a general expression for reliability, 
applied load and resistance of the structure (or system) 
of unspecified distribution are considered. ,v 

The reliability of a structural system is determined 
from the basic concept that a no failure probability exists 
when resistance (R) of the component or member is not 
exceeded by the applied load (l) or resistance > load applied. 


Let the density function for load (L) be denoted by 
f L (l) and that for resistance (R) by f a (r), Then by 
definition. 


Reliability = P (R > L) = P (R - L > 0) ... (2.2.1) 

The shaded portion in the Rig. 1,2 shows the interference 
area, which is indicative of the probability of failure. 
Figure 2.1 shows the enlarged portion of the Fig. 1.2 . 

The probability of a load value L 0 lying in a small interval 
of width dl is equal to the area of the element dL, i.e. 


dl dL 

3? (I, 0 - r <I<I b +5“) = fidfa) dli ••• (2-2-2) 

The probability that the strength R is greater than a certain 
load L Q is given by 


P (R > 1 0 ) = / f a (r).dR ... (2.2.3) 

L o 


The probability for the load value lying in the small interval 
dl and the strength R exceeding the load given by the small 
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interval dir under the assumption that the load and strength 
random variables are independent is given by 

oo 

f f R (r) dR ... (2.2.4) 

■o 

How the reliability of the element is the probability that 
the strength R is greater than the load L for all possible 
values of the load L and hence is given by 

OO oo 

Reliability = f f L (l) ( / f R (r). dR ) dL ... (2.2.5) 

— OO J, 

Reliability can also be computed on the basis that the 
load remains less than the strength. The probability that 
the strength is with in a small interval dR is 
dR dR 

P(R Q - ~2 - R - R o + F ) = f R (R 0 ) dR ••• (2.2.6) 

and the probability of the load being less than Rq is 
given by 

' P ( L < R o ' =_£ f I (l). di ... (2.2.7) 

Hence the reliability of the element (or member) for all 

the possible values of the strength R is 

cc R 

Reliability = / f R (r) ( / f L (l)dl) dR ... (2.2.8) 

The unreliability or probability of failure, is defined as 

Probability of failure = 1 - Reliability 

= P ( R < L) 

Probability of failure + Reliability = 1. 


Thus 



2.3 Reliability Computations for Normally Distributed 
Strength and load: 

Many random variables encountered in the physical 
sciences appear to be normally distributed and the normal 
distribution gives an adequate approximation to the 
distribution of many other measurable random variables. 


Thus, if a complete theory of statistical inference 
is developed based on the normal distribution alone, then 
one has in reality a system which may be employed quite 
generally, because other distributions can be transformed 
to or approximated by the normal form. 


If I and R are normal variables, then the density 
functions are 

1 

exp 


f L< ' 1 ^ S T f 2 % 


a L-L 2 

r<T) 

L 


f R (r) 


1 


S a V27t 


exp 


1 R-R 2 

2 ^ S R 


— 00 < 1 < 00 

(2.3.1) 

- 00 < R < 00 

( 2 . 3 . 2 ) 


Reliability is the probability that strength exceeds the 
load 

R - 1 > 0 

or. Y > 0 

where Y = R - 1 

Then it is well known that the random variable Y is normal: 

distributed with a mean of 
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(2.3.3) 


(2.3.4) 

and the density function is given by Fig. 2.2* 
where 

R and 1 L are the mean values of strength and load 
respectively. 

and Sjj are the standard deviations of strength and 
load respectively. 


K - li 


and standard deviation of 



S L + S R 


'Die reliability can be expressed interns of Y as 


Reliability = P ( Y > 0) 

~ 1 


/ __ 

0 S-Jf 2 tc 


exp 


1 Y-Y 2 

" 2 ( ^7 ) _ 


Y T Y 

If we let Z = ( ■ — ) tden s y dz ~ d ?« 


dY • 

(2.3.5) 


Then Y = 0 , trie lower limit of Z is given by 


0 - Y 

~“s ” 


R - 1 


• • i 


( 2 . 3 . 6 ) 


and when Y - + 00 , the upper' limit of Z -* + 00 . 

Therefore 

1 00 -7 2 /0 

e L d dZ ... (2.3.7) 


Reliabilit3 r 


■ R- Ii 


f 

R- 

i, 2 


r s 2 +s| 
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Clearly, the random variable z= (Y-Y)/S y (2.3.8) 

is the standard normal variable. 

Hence reliability can be found by merely referring to the 
normal tables (Jppendix A). 

Equation (2.3.6) which is used to find the lower 
limit of the standard normal variate Z, is commonly known as 
coupling equation. 

Equation (2.3.7) can be rewritten as 

R - 1 

Reliability = - 0 } 

= 1 - 0(2) ... (2.3.9) 

Where 0(2) is the standardized cumulative distribution 
function. 

2.4 Extended Reliability Theory: 

If R and 1 are random variables, there is a finite 
probability of failure as defined earlier 

Pf = P ( R < L) ... (2.4.1 ) 

In practice R and 1 are of course unknown, it is 
possible only to make predictions or estimates of what these 
are, or might be, using theoretical models; in particular, 
these predictions would normally be limited to the estimations 
of the means and variances of the relavent design variables.. 
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Such theoretical predictions are invariably imperfect and 
thus will contain errors. 

Assuming R and 1 as random variables, the uncertainities 
associated with the inherent randomness in R and 1 are modeled 
in equation (2.4.1 ), however, since there are also Uncertani- 
ties associated with the errors of prediction, this is not 
sufficient. To be sure, the prediction errors must be 
included in the risk model; for this purpose, the extended 
reliability model is apropos, the essence of which can be 
described as follows. 


2.4.1 Formulation of Extended Reliability Model: 

/v A 

Let R and L be the theoritical models of R and L, 

respectively; in order to compensate for any shortcomings 

in these models, apply correction factors E r and such that 

A 

R = E r R 

L = H-l 5 ... (2.4.2) 

Thus the probability of failure becomes 



= P (R < 1) 

A A 

= P (E r R < RpL ) 


(2,4.3) 


A A 

R and 1 are general ly 


random variables with mean values 


R and L , and standard deviations 35 and s£ ; the coefficients 

. = Sj 

of variation °p - -g— > u f • 
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j^ r and should also be assumed to' be random variables 
with mean values 1.0 and coefficients of variation 

. A A 

and representing the errors in R and L respectively. 

On the basis of first order theory, the mean values of 
R and 1 by virtue of equation (2.4.2) are 


S R Z R = R 
and S L L = I 

in other words, the predicted means R and h are unbiased 
estimaters of the true means S R and Sj,, where as the 
total uncertain! ties associated with the prediction of 
and randomness of R and 1 are 


2 2 2 
Vr ~ 6r +A s 

^ 2 2 

\ respectively. 


'The equation (2.4.3) can be written as 

P f = P (S < I ? ) ... (2.4.4) 

where R = R]/^,. 


R and 1 may be assumed to be independent lognormal 
variates with means R and 1 and coefficients of variation 
6 r and 6 R , whereas R r and R-j_ are also lognormal with means 
1.0 and coefficients of variation and/i^; then the 
failure probability of equation (2.4.3) is given conveniently 
by 

In R/ Jj 

(2.4.5) 


In W L 

ca. 1—0 (.^”2 — ■" fla ^f ,s ) 


* * • 



26 


where $(-) is the standard normal probability available 
widely in tabulated form. 

2.5 Observations on Reliability Theories: 

The failure probability is a function of the 

probability distributions of R and 1. So it is the 

distribution of the safety margin (R-l) or ln(R/l) that 

is important in the evaluation of the probability of 

failure than the distributions of the individual variates. 

If the risk level is large (say > 10”^) the probability 

distributions are not too important because the calculated 

failure probability would not be too different regardless 

of these distributions. However, for very small risks, 

-3 

< < 10 , the failure probability can be quite different 

depending on these distributions. In this latter case, 
knowledge of the correct distributions for R and 1 becomes 
iimoortant,. Unfortunately, with limited information, the 
correct distributions are most difficult to ascertain. 

In practice, the failure probability may only be calculated 
on the basis of assumed distributions. Consequently, at very 
low risk levels, the calculated failure probability is sme- 
v/hat s.n ambiguous measure (strictly speaking) of the 
under Ijring risk. In spite of this, however, the failure 
probability remains as the only quantitative measure of 
risk, and at the very least, a failure probability estimated 
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on the basis of a prescribed set of distributions remains 
useful as a consistent (even if only approximate) measure 
of the underlying risk. 




FIG-2-2 DENSITY FUNCTION OF RANDOM 
VARIABLE y 


CHAPTER 3 


DESIGN FORMULATION AND DESIGN CHARTS 
3 . 1 General : 

This chapter can be broadly divided into two sections 
First section deals with the parametric study of the relia- 
bility based design of RCG circular slabs. The other section 
contains formulation of design equations taking into account 
all the variabilities inherent in the material properties 
and presentation of design charts for a preassigned proba- 
bility of failure based on Indian Standard Code of Practice 

(g? ) 

for plain and reinforced concrete '. The formulation is 

kept quite general and the variables are grouped together 

to reduce their mx-iber and bring them into a non dimensional 

fornn The coefficients of variation for different design 

variables have been taken from the test results on more 

( 22 ) 

than 100 slabs compiled by MacGregor and Mir za v ' to keep 
the design as realistic as possible. Moreover to cover 
the unknown uncertainties, design equations are formulated 
taking relative ratios of the coefficient of variation of 
different design parameters in terms of a variable V. 3y 
changing the value of V a different set of coefficients of 
variation can be considered in the design, although their 
relative ratios will remain same. 
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?or calculating the external moment on the slabs 
lowerbound solutions for circular slabs are used, 3he 
design charts are presented for both simply supported 
and fixed edge conditions, Bor circular slabs with free 
and fixed edge conditions the upperbound and lowerbound 
solutions are the same. So it is more appropriate to 
consider this unique failure mode while calculating 
external moment. Only the uniformly distributed load 
is considered. 

3.2 Parametric Study of Reliability Based Design of 
■ Reinforced Concrete Circular Slabs:. 

In most of the cases the slabs are usually designed 
as under-reinforced sections. The basic reliability equation 
(2. 3.8) has been used for investigating the effect of 
different design parameters on the reliability design of 
under-reinforced concrete circular slabs. ’The parameters 
included in the present study are depth, .stress in concrete, 
stress in steel, area of steel, coefficient of variation 
and the specified probability of failure. 

'The basic aim of this study is to find the effect of 
each parameter on the design qualitatively only keeping 
all other parameters as constant. This is an investigation 
from which an educated guess to the combined effect of some 
of these parameters on the design can be thought of. In a 
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realistic problem their effect on the design considering 
the interaction (if any) should be studied to have the 
over all picture of the problem. However the present 
study has been used as a preliminary investigation before 
formulating the reliability design of concrete circular 
slabs by clubbing a number of design parameters together 
in such a way that the whole design is in non-dimensional 
form. 

Starting with an initial design, increments are 
given to each parameter turn by turn and the effect on 
area of reinforcement is found. 


3.2.1 formulation of the Problem: 

for a reinforced concrete slab the equation (2.3.8) 


can be written as 

Mr - Mr 


( 3 . 2.1 ) 


^ s m r + Sm-r 


where, 


M R = Mean expected moment of resistance of the slab 
per metre width 


M n 


Mean design moment per metre width of slab 


Sjf- and S M - =Standard deviations of M R and Mr respectively. 
X R. i Ij 

■The value of Z can be obtained from the standard 


normal distribution tables for a given value of P f 
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(probability of failure) and vice versa (Appendix A). 

Mostly slabs are designed as either balanced or 
under reinforced sections. Ming the moment of resistance 

equation from Indian Standard Code of Practice for Plain and 

, _ / Tq . AS6-1978) for under reinforced 

Reinforced Concrete (.IS. 4oo y/ ' 

section (Appendix B). 


M' 


— A._i_f ( d — K 
•R A st I y 


at *y 

f 1 t> 
ck 


) 


( 3 . 2 . 2 ) 


where 


== area of reinforcement 
d = depth of the slab 
h = metre width of the slab 
f Y = tensile stress of steel 

f , = characteristic strength of concrete 

J -ch 

K = 0.774. 

Here the partial safety factors for both f y ^ 

, k are ommitted because the variabilities in the stresses 
1 considered in the form of coefficient of variation, 
per a circular slab the lower bound solution for external 
moment with uniformly distributed load is 

p*R 2 

(% + m ) = “ 


• * • 


(3.2.3) 


where 


R 


m 

P* 


radius of slab 

negative moment, 
uniformly distributed load 



33 


let m = -^ M l 

^ere ^ is the ratio ol negative moment to positive 
moment* 

Then equation (3.2-3) can be written as 

P*R 2 . (3.2.4) 

% " 6 ( 1 +-^) 

where p = intensity of uniformly distributed load. 

It is obvious from equations (3.2.2) and (3.2-4) 

that M r and M L are the functions of a number of design 

variables Xj_. In all practical cases these variables 

are uncertain to a certain degree. These uncertaimtres 

have a considerable effect on the ejected values and 

. . ¥ an d m t . exact statistical 

the standard deviations of % and i L 

. . f m and M t can be obtained if the exact 
distribution of % aM1 L 

distribution functions of the variables X, are ,nown 

oan hope to estimate the mean and coefficients 
practice one can hop . + 

. tables. Tbe underlying probability 

of variation o _ assumed as daussian 

distribution of the random variables 

distribution, me mean and standard deviation of a fu 

oan be calculated using the mean and standard deviation 


the random variables 

Y = f (h 
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V) 


( 3 . 2 . 5 ) 
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S 


Y 


n fa f(X) 

i=i 1 ax ± 

V 


(x 1f x 


2* 



(3.2.6) 


where 

Y = function of random variables X-j^. 

Sy= standard deviation of a function of random 


variables X^ , X^. 


Bar (-) over the variable indicates the mean' of 
that variable. Here after the bar is dropped for the sake 
of convenience. 


Applying the foregoing equations (3.2.5) and (3.2.6) 
the mean and standard deviations of the functions M R and 
Mjr can be derived as follows. 


M- 


■a 


, K ^t f y \ 

A st f y ^ d b f ck 


(3.2.7) 


2 

%■ 


H 


aa R 2 a + ( ^_) 2 S | 

S *st tr V ok f 


L y 

dMp 2 2 / dM R 4 0 2 

+ <*r> s * + ( 5r“ } Sb 


fy 

2 


. (f a - ) 2 4 - 

\ X \r -P K 


f ok b 


2 

Ast 


2K Agt f y 


+ (a ,a - - . )- r. 


■ok 


’22 

y T fy 


K 4 f y >2-2 2 


f r— - } f ck Tfck 
I ck b 
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+ ( 


2 2 p 2 

< A s t 112 % 

K , 2 2 2 


f ck b 


( 3 . 2 . 8 ) 


where , V^y, Y^ c ^, V^, V^ are the coefficients 

of •variations (COY) of A g ^, f , f^, d, b respectively 

Standard deviation 

GOV = — — — 

Mean value 


The standard deviation of can be obtained from 


M t = 



6 ( 1 +^ ) 

M 1 (4 vf + V^+ T | ) 


(3.2.9) 

( 3 . 2 . 10 ) 


where, 


7 R and Yp^are the COV's of radius of slab and 
uniformly distributed load respectively. 


Since the behaviour has to be studied qualitatively 
only, the same value of Y will be assumed for all the 
variables. Thus 


7 Ast = 7 fy = 7 f ck = 7 d = 7 b 


7 R " 7 p*~ 7 * * 


(3.2.11 ) 


By virtue of equation (3.2.11), equations (3.2.8) and 
(3.2.10) reduce to 


■%R 


I 2^(f v d - ^2: ) 2 Aj 


V 

;+ 72 f y 4t d2 


2Y^K^ A g £ fy 

f ok b 

(3.2.12) 


• • • 
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2 2 
rC . tit*- 


(3.2.13) 


!:fow substituting equations ( 3.2.7 ), (3.2.9), (3.2.12) 
and ( 3 . 2 .I 3 ) in equation ( 3 . 2. 1 ), one gets 


K A i ‘ f2 

*y J ** a -r^ — * - M I 

CxC 


I 2V “ (f y d 


22 il st f y >2 


-*-ck b 


\ a 2 OTr 2 K2i 4t f y 
} A st + 2V ^^2- 


.. (3.2. H) 


+ V 2 f 2 A s 2 d 2 + 6l| V 2 


1 1/2 


Simplification of equation ( 3 . 2.1 4 ) yields the following 


fourth order polynomial in 


K 2 f,t 


st -2 k 2 
1 ck b 


*3 (10V 2 Z 2 -0 + ^ j 


3 I 2Kd 


^ok b 


P 9 N 
1 -4K7 Z ) 


+ A^ t f 2 d 2 (3V 2 Z 2 -1 ) - Mp 

™ y f ck b 


+ A st I 2 f y d M 1 


+ m| I 6Y 2 Z 2 - 1 j = 0 


A computer programe is written for solving the equation 
(3.2.15) by bisection method and Newton-Raphson method. 

The effect Of o Sher parameters is studied on area of steel 
turn by turn while keeping all other parameters constant.. 


The results are as shown in Table II and Illi Area 
of steel versus other design parameters is indicated in 
Rigs. 3* 1 &nd 3.2. 
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Starting design parameters : 
load = 6.0 kN/m 2 

Radius of slab = 5 m 
Depth of slab = 100 mm 
f y = 250 MPa 

f ck =15 MPa 

The slab considered is fixed along the edge with 
y =1.0 ; V = 0.05 ; Z = 1.0. 

TABUS II 

EFJplCT OP r f ck ' AMD *f * OM AREA OP STEEL REI MPORCEMEMT 


f ck 

(MPa) 

Area of 
Steel 

(mm 2 ) 

Percent ate f„ 

Decrease in y 

Area of Steel 

(MPa) 

Area of 
Steel 

(mm 2 ) 

Percentage 
Decrease in 
Area of Steel 

15.0 

630.82 

0.00 

250.0Q 

.630.819 

0.00 

90.0 

616.05 

2.34 

415.00 

380.011 

39«76 

25.0 

607.81 

3.64 

500.00 

315.41 

50.00 

30.0 

602. 66 

4.46 




35.0 

599.05 

. 5.04 
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TABLE III 

EFFECT OF 'V', 'Z', 'd' OB .AREA OF STEEL RSI EFORCEMERT 


Area of 
Steel 

A st 

(mm 2 ) 


Percen- 

tage 

Increase 
in A st 


Area of Per- 
Steel centa- 


A st 

(mm 2 ) 


ge 

Increa- 
se in 


d A^ Percen- 

(m “) Decrease 

in A st 


0.050 

630.82 

0.00 

1,0 

630.82 

0.00 

100 

630.82 

0.00 

0.075 

683.15 

8.29 

2,0 

740 T 23 

17.34 

110 

564.01 

10.59 

0.100 

740.20 

17.34 

3.0 

873.14 

38.41 

120 

510. 80 

19.02 

0.125 

803.09 

27,31 

4.0 

1043.42 

65.40 

130 

467 . 25 

25.93 

0. 150 

873.H 

38.41 

5.0 

1282.36 

103.29 ' 

140 

430.85 

31.70 


The following observations are made from the above study r 

i) The compressive strength of concrete has no significant 

effect on the area of steel. Using higher grade concrete 'will 
result in only nominal decrease in the area of steel. The 
curve shown in Fig. 3.1 A^ Vs f^ is almost flat. 

ii) With the increase in tensile strength of steel the 
required area of steel is decreased considerably. Almost 
50 percent of the steel area is decreased with increase in 
grade of steel from 250 MPa to 500 MPa. 

iii) The reliability parameter (i.e. probability of failure) 
is found to be more sensitive. With increasing values of Z 
the area of steel increases rapidly. 
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iv) The coefficient of variation is also a sensitive 
parameter* The increase in coefficient of variation results 
in increase of steel area. 

v) With the increase in depth of the slab the area of 
steel decreases. At greater depths the percentage decrease 
in steel area is some what less than what it is at smaller 
depths . 

3.3 Reliability Based Design of RGC Circular Slabs with 
Uncurtailed Reinf orcement : 

Based on the classical reliability approach, an 
attempt has been made in the present work to formulate the 
design concepts for under-reinforced concrete circular slabs 
in flexure only and prepare design charts for the same. To 
cover the effect of unknown uncertainities the coefficient 
of variation (GOV) is also taken as one of the variables, 
although their relative ratios have been maintained for all 
values of GOT considered on the problem. The slabs are 
assumed to be with isotropic reinforcement, which is most 
common in circular slabs. The design for shear strength is 
not considered in the present work, because in most of the 
designs for RGG slabs shear is not critical. 

In the next section the formulation of equations is 
presented. These equations will be similar to those developed 


4 


in case of parametric study. But here all the design 
variables are converted into non-dimensionalized form, 
from which an equation of fourth order polynomial in terms 
of steel area is deduced* 

3.3.1 Formulation of the Design Problem: 


Usually the slabs will be designed as under- 
reinforced sections and their moment capacity as per Indian 
Standard Code of Practice for plain and reinforced concrete 
can be obtained as follows (Appendix B): 


(32) 


M- 


R 


K A q+ f_ 

A,* f_ (d - — -2t~ y ) 
S 7 b f=k 


(3.3.1 ) 


where K = 0.774 

Balanced steel area is obtained from 




max 


T3t 7 


0.5423 f cJc db 


(3.3.2) 


wnere 


and 




max 


is maximum permissible neutral axis depth 


u,max 

d 


0.0035 


0.0055+ 

^s 


(3.3.3) 


from equations (3 . 3 « 2 ) and ( 3 . 3 *■ 3 ) 


0.0035 


54. 23 f ek 


0.0055 + fy/Bj, f y 


(3.3.4) 


p n = Percentage of balanced reinforcement. 


From the above equation one can see that the balanced area 
of steel reinforcement will depend on grade of steel and 


concrete. 



41 


P 


o 


A st 

bd 


x 100 


xu. Yxcunig equation (3.3.1) by a quantity bd 2 f cfc throughout, 
a non-dimensional expression is obtained as 


a R 


bd 2 f 


ek 




bd f ck 


bd f ck 


Taking 


Q r = 


P 


M 


•R 


bd2 f ck 

— x 100 
bd 


m 


100 fok 

Now equation (3.3.2) can be written as 


Qx> ~ pm (1 - K pm) 


... (3.3.5) 

... (3.3.6) 

•. (3.3.7) 

... (3.3.8) 

... (3.3.9) 


p ii«ab been expressed as percentage of steel. 


The external moment for uniformly distributed load 
on the circular slab is given by the following expression 


M i = 


-* p 
P 


6 ( 1 +^ ) 


.. (3.3.10) 


where P* ~ uniformly distributed load. 

Ifow converting the above equation into a non- 

2 

dimensionali zed form by dividing the equation with d 
through out* 
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# 2 
p R^ 


6(1+S))d 2 f ck 


p* R 2 

( ) ( 7 ) 


Q l = PL 2 / 6 (1 +"'■) ) 


6(1 ) 


where Q L = M L /(d 2 f clc ) 

p = pVf ck 


(3.3.H) 


(3.3.12) 


(3.3,13) 


(3.3.H) 

(3.3.15) 

(3.3.16) 


Row assuming Q R , Q R as the normal variates, the reliability- 
parameter ’ Z' from equation (3.2.1) is 


Z = 


%t - gi 

<S§ ♦ ) V2 

^R 


... (3.3.17) 


where Q r and Q r are functions of several variables. 

Q r - f(d, f ck , fyt ^- s t> 0 ) 
q l = *<p* R, a, f ok ) 

Row assuming 7 * vu . y 7~ , 7~ , 7 a , V h as the 

p K d ck £ y sst D 

respective coefficients of variation for variables p*, R, d, 
f jj., A gt , b and using the equations (3.2.5), (3.2.6), (3.3.0 
through ( 3 . 3 . 17 ) the following expressions for standard 
deviations of load and resistance can be obtained. 

2 , ac * R N 2 „2 . / 0 q Rn 2 *,2 ... ( 3 . 3 . 18 ) 


s 2 = rJSLf s 2 4- (IS) s 
% Sp ' p ' E 
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2 2 p dQ-j- 2 2 

s Q t = (rr- ) S 2 + (— y Sj 

■• • • 

(3.3.19) 

dQr) p 

(— ) = (.m (1-2K pm)) 2 

« • « 

(3.3.20) 

d Qp 2 o o 

(-—) = (p - 2 Km p 2 ) 2 
om 

• • • 

(3.3.21 ) 

o dp 2 2 dp 2 2 

s p " ( a^> St + ^ Sb + 

(-) 2 s 2 
^dd ; b d* * 


= p 2 (v| 3t + + V 2 ) 

• • • 

(3.3.22) 





2 TT 2 , 2 ,,2 
= m VP + m V- . 

7 z ck 


= m 2 (v| + V 2 ) 

I y x ok 

... (3.3.23) 

5 dp p ? , \2 q 2 

S P dp# p ^df ck f ck 

... (3.3.24) 

= P 2 (V 2 + V 2 # ) 

± ck P 

(3.3.25) 

o d Ij 2 o d Ji n o 

s? = (— ) s 2 + (— ) 2 s ^ 

p dR ' R dd ; d 

(3.3.26) 

= (~) 2 (v 2 + V 2 ) 

M 'H d ' 

(3.3.27) 


substituting equations (3.3.,22) and (3.3.23) in equation 
(3.3.18) 

Sq r = P 2 jm 2 (l_2K pm) 2 (V 2 ^ + V 2 + V 2 ) 

+m 2 lp 2 (l-2K pm) 2 (V 2 + V? ) 

|__V' ' J 7 • r ck 
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Sq = (V 2 m 2 + 4K 2 p 4 m 4 - 41 p 3 m 3 ) x 
R 


(? i + 7 k + ti 

st d d 


v£ + Y? ) 
~y I ok 


(3.3.28) 


jubs t, it using equations (3*3.25) and (3.3*27) in equation 
(3.3.19) 


q 2 


L 1 


,2 * 2 '(V?* + Vf) + 


2 2 
4P^ L 


L 36(1 p X y 35( 1+ -^)2 


O O P 

l (Vp + vr 


3,3 


r2 


(V.f* + V| ' + 4Vg + 4 .7? ) z--p 

P 1 ck K d 36( 1+ ^)2 


p 2 i 4 


(3.3.29) 


OOY's for different variables using the tests results on 

( 2 °) 

slabs by MacGregor and Sllingwood 


Vf , 

ii 

o 

• 

00 

v f 

= 0 

.107 






ck 


x y 








7 A 

st 

= 0.10 

V 

= 0 

.12 






7 d 

= 0.017 

7 b 

= 0 

.01 







= 0.01 

7 e 

= 0 

.033 

(COY 

of 

jSIs 

i.stic I 

iodulus ) 



s 






(3 

.3.30) 

■ the 

relative 

ratios 

of 

the 

COY's 

in 

(3. 

3.30) 

and 


still treating it as a variability in 7 , they may be expressed 
as 


V f = 18 I ; I* = 10,7 7 
x ck . ~y 


• y. 

» v ii.. 


10 V 


‘st 


\* =12 7 ; 7 d 


1.7 V; V R =V;7 b = V 


; ^ * 1 ) 
1 • > * • / 


• * * 
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Then from equation (3.3.30) and ( 3 . 3 . 31 ) 


4 r 


4 


= (p ^ 2 + 4K 2 p 4 m 4 - 4Kp 3 m 3 ) 

((107) 2 + Y 2 +(1.7Y) 2 +(10.7Y) 2 +(18Y) 2 ) ■ 

= 542.38 Y 2 ( pV + 4K 2 p 4 m 4 - 4Kp 3 m 3 ) ... (3.3.32) 
P 2 L 4 

36(1 +~5 ) 2 


( ( 1 2Y ) 2 + (18Y) 2 + 4Y 2 + 4(1.7Y) 2 ) 


p 2 L 4 


= 483.56 Y c 


... (3.3.33) 


36(1 + ) ; 


low substituting the equations ( 3 . 3 . 9 ), ( 3 . 3 . 13 ), ( 3 . 3 . 32 ) 
and (3.3.33) in equation ( 3 . 3.17) 

pm (1 - K pm ) - PL 2 / 6(1+'") ) 

| 542. 38Y 2 (p^m 2 +4K 2 p 4 m 4 -4Kp 3 m 3 )+ 483 . 56Y 2 P 2 L 4 /36 (1 +h'0 2 | '^ 2 
l — _~j 

... (5.3.34) 


Simp lying the above equation 

3 2 ( 542 . 38 Y 2 (p 2 m 2 + 4 K 2 p 4 m 4 - 4 Kp 3 m 3 ) + 483.56Y 2 P 2 i 4 / 36 (1 +V) 2 ) 

= pm - K p^ 2 - PL 2 /6(1+- : 0 2 
further simplication yields 

p 4 (2169.5 2 Y 2 K 2 Z 2 m 4 - K 2 m 4 ) + p 3 (2Km 3 -2169.52Y 2 Z 2 Ek 3 ) 


+ p 2 (542.38Y 2 Z 2 m 2 - 


Em 


+ P (; 


) + 


36 ( 1+^) 2 


2 PL 2 /3 ( l+"'0) 

(433. 56Y 2 Z 2 -1 ) = 0 


(3.3.35) 
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The above equation (3.3* 35 ) is a 4th order polynomial 
in p (percentage of steel) and is of the same form as that 
of equation (3. 2. 1.4)* It can be solved for p (percentage of 
steel) for different values of m, P, L, V, Z. It can be 
noted from the above equation (3.3.35) that the number of 
variables are decreased to four from seven in the case of 
equation (3.2.14). Moreover by simply changing the value 
of K from 0.776 to 0.91 design charts for 15 4-56 :1964 
can be produoed. And from changing the value of K from 
0.776 to 0.59, the above equation can be used for the design 
of slabs' according to American Gode of practice. 

The above equation (3.3.35) is formulated on the basis 

of equation from Indian Standard Gode of Practice for plain 

( 3 ? ) 

and reinforced concrete . The design is restricted to 
isotropic reinforcement and inflexure only. So the values 
of p obtained from the above equation (3.3.35) will be the 
reinforcement based on limit state of strength in flexure 
only. Design charts are presented from the same equation 
(3.3.35) for different values of f y » f ok , V, Z. The graphs 
are drawn for P versus p. p*is taken equal to a minimum of 
3 idtf/m 2 to a maximum of 12 kl/m 2 . The load range on slabs 
generally will vary between these values only. The percen- 
tage of steel is restricted to balanced percentage of 
reinf orcement for different values of fy and f ck « 


Some 
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charts are presented for both simply supported ROG slabs 
onm e for .'.ixed edge conditions. Knowing the value of 
(R/d) and load p* one can directly pick up the value of 
P (percentage of reinforcement) from the appropriate chart. 
The use of charts is illustrated in the final chapter by 
solving two examples, one each from the above mentioned 
two edge conditions at the boundary. 
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FIG-3-2 VARIATION OF STEEL AREA WITH 
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CHAPTER 4 


SOME ECONOMIC IMPLICATIONS IN REINFORCED 
CONCRETE CIRCULAR SLAB DESICN 

4.1 Introduction: 

0 

The basic aim of an optimum structural design is 
to arrive at a solution which satisfies a specified 
objective in the best possible manner. Usually, the 
specified objective functions are the minimum weight, 
minimum cost, minimum volume, maximum stiffness etc. The 
present study deals with the minimum cost design of concrete 
circular slabs with uncurtailed reinforcement for a specified 
reliability or probability of failure,. 

It is seen from the Chapter 3 that a more rational 
criterion for structural safety is the reliability. Most 
of +Vm earlier minimum cost designs consider the design 
parameters as deterministic quantities. Recently however 
some studies have been performed considering the statistical 
nature of the material strengths. 

In probability based design methods, usually a 

design is considered to be safe and adequate if the 

probability of failure of the structural member is less 

—6 1 

than a specified small quantity, say 10 . Since the 

quotation of such small probabilities may have little 
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significance, some alternative proposal for strength 
specification of structural components have also been 
made. 

It has been suggested that the design strength be 
set equal to the expected value of the strength minus a 
constant number of standard deviations of the strength. 

Such a procedure is incorporated in the present study. 

The minimum cost design of ROG circular slabs is . studied 
in some detail. The depth of the slab and the area of 
reinforcing steel are taken as the design variables. The 
variability in the material strengths, load and radius of 
the slab are considered in the form of coefficients of 
variation. It is to be noted that the approach can also be 
adopted for the optimum design of concrete beams with some 
modifications . 

The design problem is formulated as a nonlinear 
mathematical programming problem. This formulation requires 
the determination of the mean and variance of the moment of 
resistance of the slab section for one metre width. The 
constrained optimization problem is solved by sequential 
unconstrained minimization technique (SUMT). Tlumerieal 
results are presented by varying all the parameters affect- 
ing the design problem. 
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4.2 Formulation of the Optimization Problem: 

ihe design Oj_ a circular slab rs dexined by two 
variables, namely steel reinforcement A g ^. and the slab 
depth d, provided the concrete and the steel strengths 
are specified. Concrete code requirements impose limits 
on the reinforcement and slab depth. For the range of 
concrete and steel strengths available and with in the 
limits imposed on design parameters there are still a very 
large number of acceptable solutions to the design problem 
The object, then is to find that solution, which results 
in a minimum cost and still satisfies all strength and 
safety requirements. 

A. classical mathematical programming problem is 
one in which a multivariate function f (X) (where X is a 
n-dimensional vector consisting of X^, i = 1, 2,.. .n) 
is to be minimized subject to given constraints g.= (X) > 0, 
j = 1,2, .... m. The function f(X) is called the objects 
function and its choice is governed by the nature of the 
problem. In rhe present problem, this function is taken 
to represent mhe cost of the RGG circular slab per metre 
length. 

‘Taking depth of slab and steel reinforcement as 
the design variables (i = 1 , 2) in this problem the 
objective function f(X) can be expressed as 

■ ■ : .. ft 
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i(x) = c 0 .b.a + c s . A st (i+v) + Of.d ... (4.2.1) 

where Q c — cost; of ready— mix concrete per unit volume 

0 S = cost of reinforcing steel per unit weight 

= cost of thickness related form work per unit 
area 

b = metre width of the slab 
d = depth of the slab 
Ag-j. = area of positive reinforcement and 
V - ratio of negative to positive reinforcement. 

As per the reliability considerations, the slab will 

be designed for the expected moment capacity M R greater than 

or equal to the specified moment plus a constant number 

of standard deviations of the moment capacity Sm. Thus 

R 

the reliability based design criterion will be 

% > M l + G . S Mr ... (4.2.2) 

where G is a constant. 


ilow the optimization problem can be formulated as; 
Find X for a minimum value of f(X) such that 



M L - G S M A ° 


(4.2.3) 


X<«> - x. >0 (j = 1,2) ... (4.2.4) 

3 J 

X i - > 0 (3 = 1,2) ... (4.2.5) 

J J 

The upper bound on the area of steel X^^ is taken as the 


balanced steel area. 
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X, 


(u) 


W d 


where 


0.0035 f ck 

W = (0.5423 - — — T x r 2 - ) b 

( 0.0055+ gL) f y 


( 4 . 2 . 6 ) 


(4.2.7) 


The lower bound on steel area is taken as the minimum area 
of steel required according to IS Code of Practice 
(1) _ 0.12 bd 


100 


( 4 . 2*8 ) 


lower limit on depth is taken as the minimum depth required 
to control the deflection according to IS Oode^ 2 ^. 

.(1.) _ 2R 

30 


X' 


1 “ — ... (4.2.9) 

Equation ( 4 . 2 . 3 ) needs expressions for M R , Sj 


Moment of resistance (Appendix B) 

K Ao + f , r 

% = A* • V* - ~ T ~ t L ) 

J - L ck 


R 


( 4 . 2 . 10 ) 


where K 


0.776. 

P*R 2 


6(1+v) 

By virtue of equation (3.2.6), equation for Sj^ will be 

„2 


(4.2.11 ) 


3 M 


•R 


5M R 2 2 , ai Rs2 2 , r 9M R \2 s 2 

^ <t + V % ( ^ } Sf «* 


m-n 2 2 , dM R \2 c 2 

(-A) 3‘ + ~) s b 

dd d 3b 


( 4 . 2 . 12 ) 


• * • 
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(fyd 


^ Ast p v 2 A 2 „ 2 

'ck 


Ffr 1 ^ T 4 t 


, 2K Af t f v 2 p 

+ (A ^- a - T^p > 2 f y 4 

ck ,D y 


+ ( 


2 2 

K 4 * 4 


4* • *> 


)*■ f 2 p 
' *ck v f 


ck 

2 „2 


4 4 * * 2 4 ♦ ( 4 ^- P-) 2 » 2 4 


f ok ,b 


1/2 

(4.2.13) 


where V A , V f , V- , 7., V h are the coefficients of 
st 1 y 1 ck d b 

variation of area of steel, steel tensile strength, compre- 
ssible strength of concrete, depth of slab and 1 metre width 
of the slab respectively. 


If the coefficients of variation V A = V f = Zfck = V d 

st y 

= Vfc = V, i.e. same for all the variables, then, 

" 22 v 4 


>Mi 


L R 


= V 2i St ( f y d " 


f ck* b 


■y p 


. 2x2 4t 4 . p2 ,2 ,2 

+ + f y 4 

ck D 


1/2 


(4.2.14) 


Finally the constraints can be stated after nomalization 
as follows : 

+ 0 SM R 


( 


M R 


) - 1 < 0 


X, 


- 1 < 0 


(4.2.15) 

(4.2.16) 


W X, 
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(4.2.17) 

( 4 . 2 . 18 ) 

(4.2.19) 

4.3 Method of Solution of the Optimization Problem: 

The optimization problem formulated in the foregoing 
section is solved by the sequential unconstrained minimization 
technique. The constrained problem is transformed to an 
unconstrained minimization problem by appending the constraints 
to the objective function through a penalty parameter r^ . 

Thus a new objective function is defined as 

m i 

tS (1 , r k ) = £(X) - r k S __ ... (4.3.1) 

3=1 g.U) 

J 

(j = 1 , 2, .. .number of constraints), 
where 0 (X,^) is called as penalty function. The second 
verm of the equation (4.3.1) is the penalty term and becomes 
large as the vector of design variables X approaches a 
constraint. The vector of design variables corresponding 
to the minimum value of the objective function is found 
by carrying out a sequence of minimizations of the penalty 
function for a decreasing sequence of values ^(r^ <r k ). 


R 

~v 

1.2 X 1 


1i < 0 


- 1 < 0 


X, 


- 1 < 0 
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For any fixed the unconstrained minimum of the penalty 
function is found by Davidon-Pletcher - Powell method \ 
The derivatives are calculated by finite differences method. 
The initial feasible design points for the optimization 
procedure are found by a process of trial and error. 

4.4 Results: 

Several examples are taken in order to illustrate 
the effectiveness of the method described. In these 
examples, the cost coefficients are taken as follows: 

O 

C f = Rs, 60/m and for 200 mm depth 

0 S = Rs.600/kN for f = 250 MPa 

0 C = Rs* 550/m 5 for f cfc = 15 MPa., 

Reference Example : 

Results obtained are as shown in Table IY 
by taking the following initial values. 

p* = 12 kN/m 2 ,* Radius = 3 n 

Number of standard deviations 0=3 


v = 

0.05 ; 

d = 230 mm 

Ast ~ 

2 

900 mm ; 

v = 0.0 (i 

Hi 

II 

250 MPa ; 

f ok = 1 5 MPa. 


simply supported) 
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Upper bound on depth is taken as 300 mm. 

- or this example the cost of the slab is reduced from 

Rs. 278.66 to Rs. 231.04 per metre length by the optimization 

procedure. The optimum vector found is as d = 213.54 mm; 

^st ~ 537.17 mm • A g ^. should be provided in both the 
directions of the slab. 

The variation in the optimum design with respect 
to the reference design is studied by changing the various 
design parameters and the results are tabulated in Table IY. 

(i) Effect of Number of Standard Deviations 1 0 ' : 

The effect of increasing reliabilities are shown in 
Table Y. is the reliability is increased by increasing the 
number of standard deviations * G f , the expected strength of 
the slab has to increase and hence the corresponding increase 
in the cost. This is evidenced by the results of Table Y. 

(ii ) Effect of cost coefficients on the optimum design: 

With variation of individual costs of forming, steel 
and concrete the effect on the optimum design is studied. 

The results are shown in Table VI. It is observed that the 
cost of forming has the dominant effect on the optimum cost 
design. The change in cost of steel has the least effect 
on the optimum design cost out of these three costs. 
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(iii) Effect of coefficient of variation: 

The effect of coefficient of variation is also 
studied on the optimum design. It is observed that the 
cost goes up with the increase in the coefficient of 
variation value. The results are tabulated in Table VII. 


(iv) Effect of building code requirements: 

When the limitation on the depth of the slab to 
control the deflections is incorporated as a constraint 
in the minimization problem, it becomes an active constraint 
forcing d to a value larger than that required for strength 
alone and increasing the cost by as much as 30 percent. 

The limitation on the steel reinforcement is always 
satisfied at optimum and has no effect on the cost , i.e. 
it becomes an inactive constraint. 
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TABIE iv 

MINIMUM GOST OP THE REFERENCE DESIGN 

Starting values: load p* = 12.00 kN/m 2 , Radius R = 3m 

Reliability parameter 0 = 3.0, 

V ± (i=1 ,2.. .5) = 0*05 , f = 250 MPa 

f ck = 15 MPa , x{ u) = 300 mm 
Starting design vector : X 1 = 230 mm ; X 2 = 900 mm 2 

Objective function corresponding to the starting design 

vector = 278.66 




Number of Penalty 

Iterations Parameter 

X 1 

(mm) 

x 2 

(mm 2 ) 

Penalty 

Function 

MX) 

Objective 
Function 
Cost (Rs. ) 

1 

10.90630 

222.4044 

899.8530 

405.5618 

272.1902 

2 

5.45315 

221 .5650 

830.8840 

332.7741 

265.1040 

3 

2.72657 

220.2190 

600.3594 

275.8269 

240.9962 

4 

1.363288 

217.1004 

558.3658 

257.3444 

236.1283 

5 

0.68164 

213.5607 

537.1708 

237.3433 

231.2593 

6 

0.340822 

213. 5508 

537.1706 

232.7004 

231.1527 

7 

0.17041 1 

213.5401 

537.1700 

231.9265 

231.0434 

8 

0.8520x1 CT^ 21 5.540I 

537.1700 

231. 1605 

231.0434 

9 

0.852x10“ 

2 21 3.5401 

537.1700 

231.1435 

231.0434 

1 1 

0.852x10" 

*213.5401 

537.1700 

231.0434 • 

231.0434 

12 

,0.852x10“ 

'*21 3.5401 

537.1700 

231.0434 

231.0434 


TAB IE V 


EFFECT OF C ON OPTIMUM DESIGN 

Reference values : p* = 12.00 kN/m 2 , R = 3 m 

T ± (1=1,2. .. 5) = 0.05 , £ y = 250 MPa 

f ck = 15 MPa » = 500 ^ 

2 

Starting design vector : = 230 mm , = 900 mm 

Starting value of objective function = 278.66. 


Values of Objective Function Optimum X Percentage <&*?*&■} 

0 Cost (Rs.) D»«M3fta i OQ in 

— — — ™ — Objective * 

Starting Optimum XL| Xg Function 

Value Value 2 Compared to 

(mm) (mm 41 ) Ref. Design 


1 

278.66 

208.5196 

201 -.7 3 

400.88 

25.17 

2 

278.66 

222.0570 

210.08 

470.63 

20.31 

3 

278.66 

231 . 0.434 

21 3. 54 

537.17 

17.08 

4 

278.66 

238.6312 

214.38 

610.38 

14.36 

5 

278.66 

247.3270 

214.69 

701.674 

11.244 
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TABLE VI 

EFFECT OF COST COEFFICIENTS ON OPTIMUM DESIGN 

■All reference values and the starting design vector are same 
as those given in Table V, also C = 3 . 


Increase in Cost Objective 

Function 

Optimum 

Value 


Reference Design 

231. 0434 

0.00 

213.54 

537.17 

Concrete Cost 

+ 50 Percent 

271.7 094 

17.54 

200.15 

503.57 

Steel cost 

+ 50 Percent 

240.6464 

4.10 

202.46 

494.58 

Forming Cost 

+ 50 Percent 

275.3184 

19.10 

219.61 

602.86 


Percentage 

Increase 

Compared to 

Reference 

Optimum 

Design 


Optimum X 


X. (mm) X 9 (mm2) 
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TAB IE VII 

effect of coefficient of variation on optimum cost 

Reference design values are same as mentioned in Table V, 
except variation in V^ (i=1 ,2, . , . .5 ). 

Coefficient Objective Function Optimum X 


Variation 

Optimum Value 
Cost (Rs. ) 

X 1 (mm) 

Xgfriim 2 ) 

0.05 

231.0434 

213.54 

537.17 

0.10 

258.4446 

213.92 

816.30 

0.15 

329.1570 

211.37 

1617.82 



CHAPTER 5 


RESULTS AMD DISCUSSION 


5.1 General: 

The design charts presented in Chapter 3 are 
applicable to design the RCC circular slabs in flexure 
only. In presenting the charts it is assumed that the 
reinf orcement is uncurtailed and isotropic. Both fixed 
and simply supported edge conditions are considered with 
uniformly distributed load. The limit state design of 
strength in shear is ignored, because in many cases shear 
stress is not critical in the design of ROC slabs. The 
following set of values for different design variables 
were used to evaluate the results : 

Y = 0.005, 0.0075, 0.01, 0.0125 and 0.015 

Z = 1.0 , 2.0 , 3.0 , 4.0 

f y = 250 , 415 MPa ' 

f ck = 1.5 MPa 

~\) = 0.0 , 1.0 

P_x10 5 = 0.25, 0.30, 0.35, 0.40, 0.45, 0.5, 0.55 and 0.6 
R/d = 20, 25, 30, 35, 40 and 45. 
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Cherts are presented for different combinations 
of the design variables* Design charts only for few 
combinations of design variables f and f Qlc are presented 
to illustrate the design procedure qualitatively. 

Two design examples are presented for the illustration 
of the use of the design charts. 

5.2 Design Examples: ' ' • 

Example : 1 


Design of a circular slab, simply supported along 
the boundary with the following design variables and for 
different probabilities of failure given. 


Radius of the slab R = 3 m 
Dive load w-j_ = 4 kR/m 2 


f ck = 15 MPa 
f =415 MPa 
V = 0,005 


Probabilities of failure (i) 

(ii) 

(iii ) 

(iv) 


P f = 0.158655 
P f = 0.022750 
P f = 0.001350 and 
p f = 0.000032 


Pbr the given probabilities of failure, the corresponding 
Z values are obtained from .Appendix A as Z = 1,2,3 and 4. 
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Assume the slab is isotropically reinforced at the 
bottom* Assume the effective depth of the slab d = 130 mm. 
Overall depth = 130 + 20 = 150 mm 
Self weight of the slab w' d = 3.75 klf/m 2 


lotal^load on the slab p* = w d + w d 

= 7.75 kJ/m 2 

R 3000 

= — = 23.07 
d 130 

P =(--)= 5.166x10" 4 

f ck 

A R 

i7ow from figs . 3.3 and 3.4 for P = 5.166x1o“ 4 and — = 23.07 

1 

(i) Z = 1 ; p = 0.20 ; A gt = 260 mm 2 

(ii) Z = 2 ; p = 0.25 ; A st = 325 mm 2 

(iii) Z = 3 ; p = 0.31 ; A st = 403 mm 2 

(iv) Z = 4 ; p = 0.38 ; = 494 mm 2 


Thus it is seen that the percentage of steel is increasing 
with decrease in probability of failure. 


Ex ample ■ 2 : 

Design of a circular slab, fixed along the boundary 
with the following design variables and for different 
probabilities of failure given. 

A = 4 m ; w d “ 4 kl/m 2 ' 
f c k = l5MPa; f = 415 MPa 
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Y = 0.01 

Probabilities of failure (i) P f = 0.158655 

(ii) P f = 0.0 227 5 0 

(iii ) P f = 0.001350 and 

(iv) P f = 0.000032 

for the given probabilities of failure, the corresponding 
Z values are obtained from Appendix A as Z = 1,2,3 and 4« 
Assume the slab is isotropically reinforced at the bottom 
as well as at the top near the edges, let the moment 
capacity at the support is equal to that at the mid span. 
Then =1.0 

let d = 1 30 mm 
Total depth = 1 30 .+ 20 = 150 mm 


= 4 kN/m" 


; w d = 3.75 kN/m‘ 


Total weight on the slab p* =7.76 kN/V 
4000 


R 

d 


= 30.76 


130 

P* 



f ck 





Now from Pigs. 3. 18 

for 

P 

= 5. 166x1 0" 4 

and R/d = 30.76, 

and 

3.19 




2 

(i) 

Z = 1 ; 

P 

= 

0.21 ; 

Ag t = 273 mm 

Q 

(ii) 

Z = 2 ; 

P 

= 

0.305 ; 

A gi . = 396.5 mm 

(iii) 

z = 3 ; 

P 

= 

0.465 ; 

P 

Ag^ = 604.5 mm 

0 

(iv) 

z = 4 ; 

P 

= 

0.74 

= 962 mm" 



Provide the same percentage of reinforcement at the top of 
the slab up t« a length of R/3 = 1.33 m from the boundary. 

5.3 Conclusions and Recommendations: 

■4 

The following conclusions and recommendations can 
he mentioned from the study of the reliability based design 
of ROC circular slabs, 

(i) It is observed from the design formulation presented 
in Chapter 3 the area of reinforcement increased with the 
decrease in probability of failure and with increase in 
coefficients of variation. 

(ii) At high coefficient of variation it is not possible 
to design the section for less reliability, we have to design 
the section for a minimum probability of failure. 

(iii) It is recommended to decrease the coefficient of 
variation to get a low percentage of steel area even at low 
probabilities of failure. 

(iv) Even with the rich mix concretes there is only a 
slight change in percentage of reinforcement. So it is better 
to use the grade of concrete depending upon the conditions 

of exposure (environmental conditions). 
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The feasibility of designing concrete slabs for 
minimum cost with a reliability based constraint has been 
demonstrated in Chapter 4. 'The problem of specifying 
small clumbers of the order 10~ 6 for the probability of 
failure has been avoided by taking the statistical constraint 
as shown in equation (4.2.2). This constraint can be 
directly converted into a probability of failure constraint 
if the distribution functions of the strength parameters 
are known, For example, a value of C = 2 in equation (4.2.2) 
corresponds to a probability of failure of 0.0227, if the 
moment capacity Mj^ follows standard normal distribution. 

On the basis of the results summarized in Chapters 
3 and 4, it is possible to adopt some guide lines for the 
design of concrete circular slabs. It is believed that 
these guide lines will result in a design of minimum cost. 
Furthermore, the sensitivity analyses performed indicate 
that the results will remain valid for a wide range of 
changes in the cost of materials, in the light of changing 
economic conditions. 

Design Guidelines: 

(i) Choose a steel reinforcement of high strength. 

(ii) Choose a concrete strength in the medium to low 
range . 

(iii) The coefficient of variation should be as minimum 
as possible. 
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5.4 Scope for further Work: 

In the present work the curtailment of reinforcement 
is not considered, it should he considered in the reliability 
based design. The depth of the slab is assumed as uniform 
throughout, and only the isotropic reinforcement is considered. 
These two > ortho tropic reinforcement and variable depth should 
also be considered in further studies. The strength design 
in shear is also not considered in this work. It can also 
be considered in further studies where in it is likely to 
be critical. 

A. more realistic method of design includes the 
control of mix, field conditions of cuping, local product 
availability, etc. further work can be done including 
these considerations. It is hoped that the design formulation 
presented will lead a way for many more soiphisticated designs. 
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APPENDIX A 


TABLE OP STANDARD NORMAL PROBABILITIES 


z 

f ( z ) 

Reliability 

R(Z) 

Probability 
failure 'Pf 1 

0.0 

0.398942 

0.500000 

0.500000 

0.1 

0.396952 

0.539827 

0.460173 

0.2 

0.391043 

0.579260 

0.420740 

0.3 

0.381388 

0.617911 

0.382089 

0.4 

0.368270 

0.655422 

0.344578 

0.5 

0.352065 

0. 69 U 65 

0.308535 

0.6 

0.333225 

0.725747 

0.274253 

0.7 

0.312254 

0.758036 

0.241964 

0.8 

0.289692 

0.788 U 5 

0.21 1855 

0.9 

0.266085 

0.81 5940 

0.184060 

1.0 

0.241 971 

0.841345 

0.158655 

1.1 

0.217852 

0.864334 

0.135666 

1.2 

0.194186 

0.884930 

0.115070 

1.3 

0.171369 

0.903195 

0.096805 

1.4 

0.149727 

0.919243 

0.080757 

1.5 

0.129518 

0.953193 

0.066807 

1.6 

0.110921 

0.945201 

0.0547 9 9 

1.7 

0.094049 

0.955435 

0.044565 

1 .8 

0.078950 

0.964069 

0.035931 

1.9 

0.065616 

0.971284 

0.028716 

2.0 

0.053991 

0.977 250 

0.022750 

2.2 

0.035 475 

0.986097 

0.013903 

2.4 

0.02239,5 

0.991803 

0.008197 

2.6 

0.013583 

0.995339 

0.004661 

2.8 

0.007915 

0.997 495 

0.002505 

3.0 

0.004432 

0.998650 

0.001350 

4.0 

0.000134 

0.999968 

0.000032 

5.0 

0.0000015 

0.9999997 

0.0000003 





APPENDIX B 


* MOMENTS OP RESISTANCE FOR A RECTANGULAR SECTION 
(With out Partial Safety Factors) 


Without partial safety factor for characteristic 
strength of concrete the area of compressive stress 
will be 0.54238 f ck X a b. 

For equillibrium of force 

Tensile force = Compressive force 
*st f y = 0-5+23 f ok *a *> 

Here partial safety factor for f y is also not considered 


d 


A , f 
at y 


0.5423 f ok 4b 


*u 


< £v 


- limiting value 
0.0035 


K b = 


0.0055 + fy/Es 


For under-reinforced section 

% 


0.42 fy 

Mt, = A_j. . f v ( d — — r ' 

R st y 0.5423 f^v-b 


ck l 

Vfc % 


= fyd (1 - 0.774 ^ 


) 


If 


5 l. 

d 

%, limit 


> Kx 


Xtijmax ^.max 

0.5423 ' — (1-0.42 — t 


•)bd 



